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Abstract 

The effect of a random surface field, within the bimodal distribution, on the layering 
transitions in a spin-1/2 Ising thin film is investigated, using Monte Carlo simulations. It 
is found that the layering transitions depend strongly on the concentration p of the disor- 
der of the surface magnetic field, for a fixed temperature, surface and external magnetic 
fields. Indeed, the critical concentration p c (k) at which the magnetisation of each layer 
k changes the sign discontinuously, decreases for increasing the applied surface magnetic 
field, for fixed values of the temperature T and the external magnetic field H. More- 
over, the behaviour of the layer magnetisations as well as the distribution of positive and 
negative spins in each layer, are also established for specific values of H s , H, p and the 
temperature T. 
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1 Introduction 

Several authors have studied the layering transitions of magnetic Ising systems, such as Pandit 
et al. [1], Pandit and Wortis [2], Nightingale et al. [3] and Ebner et al. [4-7]. Many exper- 
imental studies have motivated these theoretical works. Indeed, several methods have been 
used to study the layering transitions in Ising magnetic films. Benyoussef and Ez-Zahraouy 
have studied these layering transitions using a real space renormalization group [8] , and trans- 
fer matrix methods [9]. Using the mean field theory, Hong [10], have found that depending 
on the values of the exchange integrals near the surface region, the film critical temperature 
may be lower, higher than, or equal to that of the bulk. Using the perturbative theory, Harris 
[11] have showed the existence of layering transitions at T = in the presence of a transverse 
magnetic field. The effect of finite size on such transitions has been studied, in a thin film 
confined between parallel plates or walls, by Nakanishi and Fisher [12] using mean field the- 
ory and by Bruno etal. [13] taking into account the capillary condensation effect. By applying 
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Monte Carlo simulations on thin Ising films with competing walls, Binder et al. [14], found that 
occurring phase transitions belong to the universality class of the two-dimensional Ising model 
and found that the transition is shifted to a temperature just below the wetting transition of 
a semi-infinite fluid [15,16]. Hanke etal. [17] showed that symmetry breaking fields give rise to 
nontrivial and long-ranged order parameter profiles for critical systems such as fluids, alloys, 
or magnets confined to wedges. Much attention has been paid to the properties of layered 
structures consisting of alternating magnetic materials. The most commonly studied magnetic 
multilayers are those of ferromagnetic transition metal such as Co or Ni. Many experiments 
have shown that the magnetization enhancement exists in multilayered films consisting of mag- 
netic layers. It was found that ferromagnetic coupling can exist between magnetic layers. From 
the theoretical point of view, great interest has been paid to spin wave excitations as well as 
critical phenomena [18-21]. The study of thin films is partly motivated by the development of 
new growth and characterisation techniques, but perhaps more so by the discovery of many 
exciting new properties, some quite unanticipated. These include, more recently, the discovery 
of enormous values of magnetoresistance in magnetic multilayers far exceeding those found in 
single layer films and the discovery of oscillatory interlayer coupling in transition metal mul- 
tilayers. These experimental studies have motivated much theoretical work. However these 
developments are applied to a large extent powered by "materials engineering" and the ability 
to control and understand the growth of thin layers only a few atoms thick. However, Many 
experimental studies have shown that the magnetisation enhancement exists in multilayered 
films consisting of magnetic layers. On the other hand, we have shown in one of our earlier 
works [22] the existence of layering transitions under the effect of a variable surface coupling. 
Moreover, for an Ising film with a wedge, we found in Ref. [23] the intra-layering transitions 
under the geometry effect consisting on the existence a wedge. When the film is subject to a 
random transverse magnetic field [24], we found the layer-by-layer transitions when increasing 
the concentration p above a critical value p c (k) for each layer k. The aim of this work is to 
study the effect a random surface magnetic field on layering transitions of an Ising thin film, 
using Monte Carlo simulations. The paper is organised as follows. In section 2, we describe the 
model and the method used: Monte Carlo (MC) simulations. In section 3 we present results 
and discussions. 

2 Model and Monte Carlo simulations 

We are studying a magnetic thin film formed with N — 4 layers coupled ferromagnetically. 
Each layer is a square of dimension N x x N y — 64 x 64 spins. N x and N y stand for the number 
of spins in the x and y directions, respectively. 

A preliminary study showed that, when performing Monte Carlo simulations under the Metropo- 
lis algorithm, we note that the relevant calculated quantities did not change appreciably for 
small film thickness: N = 3, N — 4, N = 5 and N = 8 layers; and when varying the number 
of spins of each layer from N x = N y = 32 to 128. Taking into account the above considera- 
tions, in all the following we will give numerical results for a film thickness N = 4 layers, and 
N x — N y — 64 spins for each layer. 

A surface magnetic field H s is acting only on spins of the surface k — 1. 
The Hamiltonian governing the system is given by 

H = -J £ SiSj-Y^HiSt (1) 

<i,j> i 
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where, Si(l = = — 1,+1 are the spin variables and the interaction between different spins 
is assumed to be constant. The total magnetic field Hi applied on each site % of the layer k, is 
distributed according to: 

TT _ j H + H si for all sites i e surface k=l , . 

Hi ~\ H for all sites i e layers k=2,3,4. ^ ' 

H is an external magnetic field, and H si is the surface magnetic field distributed according to 
a bimodal law: 

V{H si ) = P 5(H S - H si ) + (1 - p)5(H si ), (3) 
so that the surface magnetic field is: 

jj \ H s with the probability p , . 

1 with the probability 1-p 

where < p < 1. For p = 1, H s is acting on all sites of the surface; whereas p = corresponds 
to a situation with no surface magnetic field: H s — 0. 



3 Results and discussion 

We present in Fig. 1 the ground state phase diagram in the plane (H, H s ) for a film with 
N = 4 layers, with different line transitions and configurations. For H s > H sc = N/ (N — 1), the 
system transits from the configuration N to the configuration 10 N ~ l at the line H — 1 — H s ; 
and the configuration lO^" 1 transits to 1 N at a constant magnetic field H = —1/(N — 1). 
Whereas the transition from the configuration N to 1^ is located at H = —H s /N, provided 
that the surface magnetic field is positive and kept less than its critical value: H s < H sc . It 
is worth to note that all these configurations are located at negative values of the external 
magnetic field H. In this work, the temperature fluctuations will not be taken into account. 
Hence, the temperature will be fixed at the value T = 3.75 in all the following. In order 
to outline the effect of the surface magnetic field concentration, we plot in Fig. 2 the phase 
diagram of the studied system in the plane (H,p), for H s = 1.5. We show the existence of 
two critical probabilities: pf and p£. The former is the layering concentration corresponding 
to the transition of the surface k — 1, and located at about R3 0.40. The surface line transition 
is terminated by the endpoint C\ located at the critical concentration p c (l) ~ 0.60. The latter 
is the film transition located at approximately R3 0.55 as it is shown in the inset of Fig. 2. 
It is worth to note that the internal layers of the film transit at different lines with different 
endpoints. Indeed, the second layer k = 2 transits at line terminated by the endpoint C 2 
located at « 0.94, the third layer corresponds to the line with the endpoint C3 located at 
~ 0.92, whereas the last layer k — 4 at the line with the endpoint C4 corresponding to the 
critical concentration ps 0.88. 

These critical probabilities depend strongly on the surface magnetic field values, at a fixed 
external magnetic field as it illustrated by Fig. 3. When increasing H S) the critical concentration 
of the surface p c (X) as wen as f° r internal layers k — 2, 3, 4, p c (k) decreases when increasing H s 
more and more. In all cases the critical concentration of the surface p c (l) is always found to 
be smaller than p c (k),k = 2,3,4 of the other layers. This is due to the fact that the surface 
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transits, in all cases, before the other layers. This is outlined in Fig. 4, for a small concentration 
p = 0.30 of H s for a fixed surface magnetic field H s = 1.5. To complete this study, we give in 
Figs. 5a, 5b the distribution of spins of the surface before and after the transition shown in Fig. 
4, for H s = 1.5 and p = 0.30. Fig. 5a shows the map of negative spins for H = —0.20 before 
the transition, whereas Fig. 56 corresponds to positive spins for H = —0.05 after the transition 
located at an external magnetic field H ~ —0.08. This magnetic field value depends not only 
on the concentration p but also on the surface magnetic field H s ; and decreases when increasing 
the concentration p. Indeed, the random distribution of H s is responsible on the formation of 
small and randomly distributed islands of negative and positive spins present in the surface 
and also in the deeper layers k = 2,3,4. When p — > 1.0 the negative spins, or positive ones, 
are assembled on lager islands in a given layer and the randomness of H s is lost. 

4 Conclusion 

We have studied the effect of a random surface magnetic field, within the bimodal distribution, 
on the layering transitions of a spin-1/2 Ising thin film formed with N = 4 layers, using Monte 
Carlo simulations. We showed that the layering transitions depend strongly on the concentra- 
tion p of the surface magnetic field, for a fixed surface and external magnetic fields. Moreover, 
the critical concentration p c {k) of each layer k decreases with increasing the applied surface 
magnetic field, for a fixed value of the external magnetic field. In order to complete this study, 
we have investigated the magnetisation behaviour for each layer as a function of the external 
magnetic field for fixed values of T, H s , and the concentration p. We have also shown the 
existence of islands of negative spins surrounded by aggregation of positive spins. 
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Figure Captions 

Figure 1.: 

The ground state phase diagram in the plane (H, H s ) for a film with N — 4 layers, with different 
line transitions and configurations. 

Figure 2: 

Phase diagram in the plane (H,p) for a fixed temperature T = 3.75 and a surface magnetic 
field value H s = 1.5, showing the layering transition probability of the surface pf as well as 
the endpoints k — 1, 2, 3, 4. Inset, the layering transition probability of the film p^ and the 
layering probability of the remaining layers. 

Figure 3: 

Critical probabilities p c {k) for each layer, as a function of the surface magnetic field H s for an 
external magnetic field H = —0.20. 

Figure 4: 

Magnetisation profiles as a function of H for different layers k = 1,2,3,4 of the film, for 
H s = 1.5 and a concentration value p = 0.30. The number accompanying each curve is the 
layer order. 

Figure 5: 

Distribution of spins of the surface k = 1 before and after the transition for a film with N = A 
layers, T = 3.75, H s = 1.5 and p = 0.30: (a) negative spins for H = —0.20 before the transition, 
(b) positive spins for H = —0.05 after the transition. 
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Fig. 2 




Fig. 3 




